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Abstract 

The decay rate of a (neutral) Higgs boson {H) going into a longitudinal VT-boson 
pair is calculated by including one loop radiative corrections in a two-Higgs doublet 
model. It is assumed that the Higgs boson H is much heavier than the W boson 
and a full use has been made of the equivalence theorem. A possibility is explored 
extensively that the precise measurement of the decay rate of H could be useful as 
a probe of the other neutral as well as charged Higgs bosons through potentially 
important radiative corrections. It is pointed out that, for some choice of the pa- 
rameters, the radiative corrections are sizable and that there exist a possibility of 
extracting useful information on the other Higgs boson masses. 
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1. Introduction 



In the last two decades there have been a lot of efforts towards discovery of Higgs 
bosons both theoretically and experimentally. (Sec Rcf. [1] for a review.) There 
is potentially a hope on the theoretical side that, even before direct detections of 
the Higgs bosons, one might be able to get indirect signatures of the existence of 
the Higgs particles: i.e., the study of radiative corrections and precision tests of 
electro-weak processes have been often motivated to see low-energy manifestations 
of unknown very heavy particles. 

Veltman [2] has once argued that in the minimal standard model with a sin- 
gle Higgs boson doublet, the Higgs boson mass dependence of radiative corrections 
would be at most logarithmic on the one-loop level and therefore the internal loop 
effects of the Higgs boson are rather elusive. Einhorn and Wudka [3] studied ra- 
diative corrections to gauge boson propagators in a general way and strengthened 
the Veltman's claim. They realized that, due to the custodial SU (2) symmetry in a 
single Higgs scheme [4] , the Higgs mass dependence is reduced to a logarithmic one 
to all orders in perturbation thoery. This fact is sometimes referred to as Veltman's 
screening theorem. 

The analysis of Einhorn and Wudka tells us that radiative corrections other than 
those to vector-boson propagators may have power-type terms w.r.t. Higgs masses. 
It is also clear that, in a model where the custodial symmetry is not respected, 
there is no compelling reason a priori that Higgs bosons do not produce power-like 
corrections to gauge boson propagators. Two-Higgs doublet models belong to such 
examples and the calculation performed by Toussaint [5] shows in fact that the 
screening theorem does not apply to the two-doublet model. This fact prompted 
several authors [6-7] to study radiative corrections to the p -parameter and to the 
muon decay constant in a two-doublet model. Analyses along the line of Peskin and 
Takeuchi [8] are also available [9-10]. All of these analyses are, however, restricted 
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to the oblique-type corrections. It is now apparent that other processes should be 
investigated in the two-doublet model as much as possible from the view point of 
the internal loop effects of Higgs boson masses. The two-doublet model could allow 
for large radiative corrections and would provide us with useful information. 

On the experimental side, it has been considered that longitudinally polarized 
gauge boson scatterings, i.e., W^^W^ — > W^Wj^ , W^Wj^ — > ZlZl, are suitable 
processes to discover Higgs bosons in future hadron colliders [11]. An elaborate 
calculations of these reactions have been carried out in the minimal standard model 
(with a single Higgs doublet) including higher order radiative corrections [12-14]. 
The decay process i/(Higgs) ^l^l constitutes sub-diagrams of the gauge boson 
scatterings and the decay width has also been evaluated in the minimal standard 
model. (For experimental aspects of measuring the decay width, see Ref. [15].) 

Considering the importance of the l^VF-scattering and the decay process, we 
have decided in the present paper to undertake the calculation of the Higgs decay 
width into a W^^Wj^ pair in the two-Higgs doublet model including loop corrections. 
Our theoretical motivation is to scrutinize the dependence of the width on various 
Higgs boson masses and explore the possibihty of getting indirect signature of Higgs 
bosons other than H, which are yet to be discovered. There are two CP-even neutral 
Higgs bosons {H, h ) and a CP-odd one {A) together with a charged one (G^) in the 
two Higgs doublet model. In the present paper, we will assume that the iJ-boson is 
the lightest among these scalar particles but is still much heavier than the VF-boson. 
We will study the internal loop effects due to the other heavier scalar particles, h, 
A, and to the decay process H — > Wi^W£ . 

There have been several attempts to get information on the masses of these 
bosons. According to the tree unitarity analysis [16] of the type of Lee, Quigg 
and Thacker [17], these masses are bounded from above as mn < SOOGeV, nih < 
TlOGeV, mo < 870GeV, < 1200GeV. This is a criterion of the vahdity of 
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perturbation theory. Similar bounds are also obtained on the basis of triviality 
arguments [18]. Suppose that the lightest Higgs boson H has been discoverd in a 
future collider. The reasonable question to be raised thereby will be whether one 
could glimpse into the existence of another Higgs boson in the mass range mentioned 
above by looking at the decay width of the discoverd H. In such a situation the 
calculation in the present paper will become very useful. 

This paper is organized as follows. We describe the two-Higgs model in Sec. 2 
and discuss the outhne of our calculation. We then proceed in Sec. 3 to calculate the 
tadpole diagrams and various two-point functions. In Sec. 4 we give the calculation 
of the Higgs decay vertex function. Numerical analyses of the decay width formula 
are given in Sec. 5. Sec. 6 is devoted to summary and discussions. 



2. The Two-Higgs Doublet Model 

Let us begin with the Higgs potential, consisiting of two Higgs doublets, $i and 
$2 with Y — 1. The most general SU (2)^ x [/ invariant Higgs potential becomes 
[19] 

$2) = -1^1 1 $1 r -1^1 1 $2 r -(/^?2*i*2 + i^u^Ui) 
+Ai I $1 1^ +A2 1 $2 1^ +A3 1 $1 n $2 r 

+A4(Re$J$2)' + A5(Im$I$2)'- (1) 

It is imperative to avoid the flavor changing neutral current and henceforth the 
discrete symmetry under $2 ~^ —'^'2 [20] has been assumed except for the soft 
breaking term (/x^2^i^2 + A*i2^2^i)- This soft breaking term may be important 
in connection with a new source of CP- violation for baryon genesis [21]. Inclusion 
of this term, however, will render our calculation clumsy to some extent, because 
of the mixing of CP-even and odd states [22]. Only for this this reason, we set 
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A^i2 = /^i2 = throughout the present paper. Inclusion of mixing between CP-even 
and odd states will be discussed in our future publications. 

In order to see the particle contents and mass spectrum, we rewrite the Higgs 
potential in terms of the parametrization of the Higgs doublets 



(2) 



V TS^^^ + hi + iZi) J 

where the vacuum expectation values vi and V2 triggering the spontaneous break 
down, are assumed to be positive without spoiling generality. 

The mass term in Eq. (1) may be diagonalized by introducing two kinds of 
mixing angles, a and (3 in the following way 




cos a — sm a 
sin a cos a 

cos P — sin /9 
sin 13 cos P 

cos /? — sin /? 
sin P cos P J 

The mixing angles are determined in fact by the following relations. 



(3) 



(4) 



(5) 



tan 2a — 



(As + K)viV2 



tan p — — 



(6) 



(7r/2 > a > — 7r/2,7r/2 > /3 > 0). The Nambu-Goldstone bosons to be absorbed 
into the longitudinal part of and Z are denoted by w and z respectively. Other 
fields h, and H are neutral while G is a charged one. The five coupling constants 
in Eq. (1) are expressed by the masses of these scalar particles together with the 
mixing angles [23] 

1 



Ai 



2v'^ cos^ p 



^m^cos a + m^sm aj. 



(7) 



''^2 = r, o^- 2 o i'^h sin^ a + m]j cos^ a) , (8) 
Zv^ sm p 

sin2Q; , 9 , , 2m^ 

^3 = 2 • o^ K - + 9 
sm 2p 



2to; 



A4 = ^, (10) 

2 

A5 = ^(m^-m^), (11) 



where v = yvf + V2 ^ 246GeV. Incidentally, masses and in Eq. (1) are given 
on the tree level by the vacuum stability condition 

IJ'l = ^ivl + \{^3 + )^4)vl 1^1 = X2vl + ^{X3 + X^)vl (12) 

Eqs. (7)-(12) define the change of our initial seven parameters (Aj(i = 1 — 5), /i.^, /i|) 
into the set {rrih, mu, mc, rriA, a, P, v). 

We are now interested in the decay process H + W£. We will use 

the equivalence theorem [23, 17] which states that if 171^ :$> M^r, we are allowed 
to replace the external longitudinal gauge bosons by the corresponding Nambu- 
Goldstone bosons. In our case, therefore, the evaluation of the process H w'^w" 
is our central concern. The tree-level interaction term dictating this process is 
extracted from the potential (1) as 

^Hww = — — sin(Q; — P)Hw^w. (13) 

One might suspect that the counter terms corresponding to this part of the potential 
are just obtained by varying the coefficients in (13), that is, by putting m% 
mjj — Srn^, v ^ v — 6v, a ^ a — 6a, and /3 — > /9 — (5/3, together with the wave 
function renormalization. This, however, does not give us all the counter terms. 
There are other terms coming from the state mixing. Changes of mixing angles 5a 
and Sp induce a mixture of the paired fields, H '(-^ h, and w G, respectively. This 
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indicates that some of the counterterms are obtained from the hw^w, HGw^ , and 
w vertices 



^hww — cosfa — I3)hvivw, (14) 

V 

2 2 

J^HGw ^ - ""^ ~ cos(a - (3){Hw^G + HwG^), (15) 



by replacing h —SaH and G 5j3w. 

After all we conclude that the counter terms for the w vertex take the 
following form 

^YYt^j Su j 

S^Hww — { 2 1 1" ( V ~ 1) + i^w — ijf^Hww 

,2 ^2 



^ba— cos(q; - ^)Ew I w 

V 
V 



(16) 



Thus to evaluate the radiative corrections to this process, we have to know Smjj, 
6v, da, SP, together with the wave function renormalization constants Zyj and Zh- 
The calculation of 5v/v is fascilitated by considering the renormalization of the 
1^-boson mass M^, i.e. 6v/v — 5M^/2M^. This has been computed previously 
by Toussaint [5] and we just quote his results, 

SM^ 1 A . r, „ , „ . m?rm\ , 

-r# = J - 2m^ + m| + m^ + m^) + „^ o — T 

(47r)2i;2 '-2^ ^ ^ ^' m\-m% m\ 

+ cos^(a - /?) In ^ + s\u\a - P) In 

(17) 



Note that the logarithmic terms are all negative. Since m\ > mjj, SM^ is minimized 
at sin^(Q; — P) — 1. 
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3. Tadpole Diagrams and Self- Energies 

Before launching into the computation of the radiative corrections to w 
vertex, we prepare the counter terms (16) by evaluating (5m^, 5a, Sf3, Z^, and Zh- 
First of all, the condition of the stability of the vacuum (12) must be replaced on 
the loop- level by including the tadpole diagrams in Fig. 1. The vacuum expectation 
values vi and V2 are determined in our case by 

9 , 9 1/x X \ 9 sin a ^ cos a 
fJ'i = ^ivl + 7;(A3 + X^)vl + Th Th , (18) 

9 X 9 l/> X \ 9 ^ COS a ^ sina 
= ^2vl + 7^ A3 + \a)vI - Th Tn . 19 

Here tadpole contributions of Fig. 1 are given by 

^ 3m?T , sin^ a cos^ a x ^ / o ^ 

Th = ^ 

2v cos p sm p 

r m?r , sin a sin^ /3 cosacos^/?, mi . , ^x-, ^ / 9n 

+{^ ^ + ^^"^ - }/i 

2v cosp smp V 



, m?T , sin a sin^ /? cos a cos^ /? , 2m^ . , -, ^ / 9 s 

+{— ^ + ^ sm a - /3 }/i , 

V cosp smp V 



3ml cos^ a , sin^ a 2n 

1 , 2 ,1 2Nsin2a x / 2 a 

— imH- H — mu) cosia — p) fi(mrr) 

i;^ ^ 2 sin 2/3 ^ ^^''^ 

r m? , cos a sin^ /3 sinacos^/9, m4 , ^x-, ^ / 9n 

^ + + " - /5 

2t; cosp smp v 

, mi , cos a sin^ /3 sin a cos^ /3 , 2m^ , -, ^ / 9 x 

-{— TT^ + ^-5-^ + COS a - /3 }/i . 

V cosp smp V 



The definition of the functin f\{w?^ is given in Appendix A. 



(20) 



(21) 
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Let us next turn to the self-energy diagrams of w. The two-point function de- 
picted in Fig. 2 becomes 2x2 matrix due to the mixing with charged Higgs boson 
G. The two-point functions are decomposed into three terms according to the three 
diagrams in Figs 2(a), 2(b) and 2(c), respectively, 

^Up') = n« + u^^i + (22) 
u^aip') = n« + nS + ^Sl (23) 

Internal particles in Figs. 2(a) and 2(c) are tabulated in Table 1. Our straight- 
forward calculation shows that the sum of the three terms are summarized in the 
following form 

n^^(p2) = n^^(p2) - n^^(o), (24) 

n^G(p') = tl^.G{p^)-tln,G{0). (25) 



The explicit forms of Tl^wip^) and Ily^G{p'^) are given in Appendix B. Eqs. (24) and 
(25) show the existence of a massless pole corresponding to the Nambu-Goldstone 
boson. For a single Higgs case it has been known that the calculation of the self- 
energy may be made easier by a manipulation concceived by Taylor [25] . In our case, 
however, we did not develop similar tricks and just simply added up all Feynman 
diagrams to reach Eqs. (24) and (25). 

The wave function renormalization of the Nambu-Goldstone Boson is finite and 
turns out to be, for the on-shell renormalization, 

Zw = l-774-^{J(2m^ + ^l + < + ^A) + ^?^^ln^ 
[inj^v^ 2 — rriQ m\ 

+ cos^ {a - (5) „^ ^„ In — ^ + sin^ a - (5) „^ % In ^} . 

mjj — rriQ mjj mf^ — itiq mf^ 

(26) 
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The renormalization of the mixing angle (3 is pecuhar if we set the renormahza- 
tion condition at = 0. (A different renormalization scheme has been used in Ref. 
[26].) The vanishing of the off-diagonal part Hyjoip^) at = tells us that there is 
no corrections to the mixing angle (5. In other words the relation tan/3 = ^2/^1 is 
preserved to higher orders 

513 = -^n^G(O) = 0. (27) 
niQ 

The last term in Eq. (16) thus gives no contribution to our calculation. 

The other renormalized quantities w.r.t. the H field are given similarly in terms 
of the two-point functions, 

Zh = 1 + YI!hh{0. (28) 
6m]j = Re(n^H(m^)), (29) 
5a = ^ Re(n,^(m^)). (30) 

It is therefore an urgent task to evaluate HuHip^) and YluHip^)- The relevant 
Feynman diagrams are again those in Fig. 2 whose internal lines arc explained in 
Table 2. All of the results for these self-energies are given in Appendices C and D. 

The mixing phenomena between H and h may be understood in the following 
way. The H quanta contain the h component due to the mixing YI}ih{p^)- Therefore 
H quanta are able to transform into h and then couple to w^w pair with the strength 
determined by (14). Such mixing phenomena are already included in (16) through 
the the 5a term (times m^cos(Q; — j3)). 



4. Loop Corrections to the Higgs Vertex 
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We are now in a position to evaluate the loop effects to the w vertex as 
shown in Fig. 3. Hereafter the Nambu-Goldstone bosons are put on the mass shell 
(Pi — pI — but considering future apphcations to Wi^W£ — > W^Wj^ we will 
keep the momentum p carried by H boson off-shell and will eventually put on the 
mass-shell {p^ = rn\j) to get the decay rate. 

There are a lot of Feynman diagrams to be calculated. We would like to classify 
the diagrams into three groups accoding to the three types in Fig. 4, that is, 

r(p2)^rw + r(2) + r(3). (31) 

These three terms come from Figs. 4(a), 4(b) and 4(c), respectively. The species of 
the internal lines in these diagrams are summarized in Table 3. 

The sum of all the diagrams in Fig. 4 is rather lengthy, but is given here for our 
use in numerical calculations. As to the contribution of Fig. 4(a), we obtain 

.6 



PW = -^sin=^(a-/3)^(p^0,0,ml,) 

4 2 

mtmi 



sm(a — /3) cos {a — P)g{p , 0, 0, m^) 



1 r 9 / cos a cos^ 3 sin a sin^ 6. ^ ■? . , 

+ — W ^ n — - sm a - (3)} 

sm p cos p 

x{{ml - rriaf sm^{a - (5)g{p^ , ml, m^, ml) 
+(mjj - mlif cos^{a - (3)g{p^, m^, ml, mjj) 



+ {m\ - m^fgi^p^, m%, m%, m\)} 

1 , 2 / cos^ /3 cos q; sin^/3sinQ; ^ / r>\-. 

+ 2^^"^-^^h^^ ^^^-2"^^^^"("-^)> 

^, / 2 2\2 I 2 2 2 2 

X (m^ - m(j) g{p , m^, m^, m^ 

IK! 

^3 

t^h{t^h - ^g) , _ cos^(^a - I3)g{p^, 0, m^, m|) 



^ mliml ml){ml ml) ^.^^^ _ ^^^^^^ _ ^^^^^2^ ^2^^ ^2^ 



^3 



+^ sm' a - (3) — — -^)g{p\ mjj, mj,, 

sm p cos p 
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t;^ sin 2/3' 

^K-m|)(ml,- m^)(2ml, + mD _ _ ^fj^ 

sin 2/5 

/ 2 2 2 2 \ 

miimif + 2m'i) , ^, ^, sin2a , , , 

^ _ ^) ^Qg2(^ _ ^)_^^(p2^ ^2 ^ ^2 ^ 



(32) 

The Feynman integral is expressed in terms of the function g{p^i ml, m|, rri^ defined 
in Appendix A. As is shown there, this function is a combionation of the so-called 
Spcncc function. (Sec Ref. [27] for a concise review of the Spence function.) It is 
almost straightforward to evaluate F^^^ numerically on computer. 

The diagrams in Fig. 4(b) are expressed by the function /2(p^, m\,m\) in defined 
in Eq. (38) and provide us with 

p(2) ^ ^^{m?cos2(Q;-/3) +m^sin^(Q;-/3)}sin(Q;-/3)/2(p^0,0) 

1 , 2 /Cosacos^/? sinasin^/?, „ 2 • / n\^ 

:;{mfj[ — ) — ZmrSima — p)\ 

v^^ sin/3 cos/3 ' g \ t-)S 

x{2m^sin^(Q; - /3) + 2m|^cos^(Q; - /3) 

I 2q; ,2 2\i 21^/2 2 2\ 

sh^^"^'^ - "^//) + ^^A}f2{P , m^., mc) 

1 r 9 / COS a cos^ /3 sin a sin^ /3 , ^ 9 . . -, 
--{< ^ - 2mi sm a - /3 } 

v"^ sm /) cos /) 

1 1 

x{-m^sin^(Q; — /3) + -mf^ cos^(q; — (3) 

I 1 2a 2 2\| 2ta'/2 2 2\ 

2 shi^^"^'' - "^//) + ^gIMp , "^a) 

~~3 ("^^ ~ - "^a) - /3) cos^ia - /3)/2(p^ m^, 0) 

~~3i'^h - ^H)i^H - ^g) sin(a - /3) cos2(a - /3)/2(p^ m^, 0) 
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2^3 



+ 2mG) cos\a - P) 



sin /3 cos /3 



sin 2/5 
)f2{p\rnH,mH 



2 , COS a sm ct 2 2 2 



1 , , 2 2 ^ sin 2a 



- — {(m^-m^)-T 



sin 2/5 



2ml} {ml + 2ml) 



I rt\ 2/ , sin 2q; „ , 2 2 2\ 
X sin(Q; - /3) cos (a - ^)^^h{v , tuh, m^) 

Sill 

l r// 2 2 NSin2o; . 2\ • 2/ /o\ I 2- 



2v^ 



sin 2/3 



+ 2mc) sin^(Q; - /?) + m^} 



x(m^ + 2ml) sin(a - /^)|^/2(p^ ml, ml). 



(33) 



Finally we give the results of Fig. 4(c) which are multiplied by 2, because we 
have to include symmetric diagrams as well 

r(3) ^ -'^{(rnl-ml)'^ + 2ml}sm{a-(3)cos\a-(3)MO,mlO) 



sin 2/3 



2m 



, sin 2a 



+ ^7^{{ml^— + 2ml) sin(a - l3) cos'(a - l3) 



-\-m 



2 - sin^ a cos^ a 



^ cos 13 sin /3 
2{ml-ml) sin 2a ^ ^ 



-{ 



sin 2p 



)sm\a-(3)}f2{0,rnl,0) 
ml sm'^{a - P) + mjj cos^{a - (5)) 



-ml cos(2q; — 2(3)} sm.{a — (3)f2{0, ml, ml 



TTijj — mQ^. 2 sin 2a 



sin 2/3 



+ 2ml) sin(2a - 2/3) 



2 / COS a sm a , , ^, , , ^ 2 2 \ 

-\-2mfj[— — x) cos(a — fj)} cos(a — /))/2(0, m^, m^) 



+2 



sin/3 cos/3 
(mi - m^)2 



sin(a - /3)/2(0, mG,m^). 



(34) 



Putting all the calculations together we arrive at the decay width formula 

1 1 



4M^ 
mjj 



\M(p'^ml)\\ 



(35) 



IGtt mn \ 

Here the invariant amplitude is given through the one-loop order by the following 
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sum 

M{p^) = r(p2) + J cos(q; - P)IlhH{m]j) - ^ sin(Q; - (5)IlHH{m]j) 

+{- + \WHH{m]j) + Z^}^ sm(a - (3). (36) 

Note that, while r(p^), n/ii^(m^), and TlHH{ni\) are all divergent, the combi- 
nation of these with the weight in Eq. (36) is finite. Although this is guaranteed by 
the renormalizability, we have confirmed the finiteness explicitly by hand. This is a 
non-trivial check of our calculation. 

Suppose that ma is much larger than all the other masses. It then turns out 
that r(p^) will be dominated ( on the dimensional account) by the terms of the 
form ttiq/v^ with possible logarithmic corrections. The same things also happen 
to the two-point functions IlHHip^) and UhHip^), which will be principally given a 
form proportional to tuq/v^. Thus the decay width is potentially very sensitive to 
the choice of mo, provided that mo is large. This situation is in contrast to the 
standard model , where the effect of the scalar boson is always veiled. Of course the 
sensitivity to ma depends upon the values of a and (3, and it is interesting to see 
the behavior of the decay width for various choices of a and /3. 

5. Numerical Analysis of the Decay Width Formula 

Now let us analyse the decay width formulae (35) and (36) by putting numbers 
into the parameters. In principle we should keep an open mind to look at every 
corner of the parameter space. There have been a lot of efforts to constrain the 
parameters in the two-Higgs doublet model from phenomenological analyses, and we 
will give due considerations for an ccconomical purpose to reduce the large parameter 
space. 

Constarints on the mass of the charged Higgs boson (m^) and the mixing angle 
P has been discussed [28] by considering the low-energy data relating to neutral 
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meson mixing ( — K^, — D^, — ) and CP- violation parameters. It 
has been argued that the low-energy data exclude small tan/3 and light ma (say, 
tan/3 < 0.3-0.4, mo < 200 GeV). In the following we will examine the decay width 
formula by setting tan P — 2 and 10 as tentative values. These values are within 
a pcrturvativc region w.r.t. the tbG coupling. Wc will also vary the charged Higgs 
boson mass as 400 GeV < mc < 900 GeV. For the values in this region it has been 
known [6] that the p parameter and the radiative corrections to the muon decay 
constant do not contradict the present experimental values. Since the H boson is 
assumed to be the lightest, we fix mn — 300 GeV and also we take mh — 400 GeV. 
These values satisfy the tree unitarity constraint as mentioned in Introduction. 

As to the choice of the angle a, there does not seem to be a thorough phenomeno- 
logical analysis comparable to the case of (3. Alternatively, the angle a may be fixed 
in a theoretically oriented way. In the minimal supersymmetric model for example, 
the parameter a is not independent but expressed in terms of the Higgs boson and 
gauge boson masses. If one takes the charged Higgs boson mass as an independent 
parameter and assumes an extremely large value for it, then sin^(Q; — /?) approaches 
to unity. In this case the strength of the w vertex (13) becomes maximal. 
In the following we will examine two extreme cases, namely, sin^(Q; — (3) — 1 and 
sm'^{a — P) = 0. Note that, in the latter case, the i/wt w vertex vanishes on the tree 
level, and the decay process proceeds only due to loop effects. 

To sum up we will evaluate the decay width (35) as a function of ma for the 
following cases. 

Case 1: tan a = tan /3 = 2, m// = 300 GeV, mn = 400 GeV, 
Case 2: tan a = tan /? = 10, mH = 300 GeV, = 400 GeV, 
Case 3: tan /3 = 2, sin'^{a - /3) = 1, mjj = 300 GeV, mu = 400 GeV, 
Case 4: tan f3 = 10, sin'^{a - /3) = 1, = 300 GeV, = 400 GeV, 
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As to the mass m^, we will set tentatively = 350 GeV, niA = 700 GeV, and 
niA — 1000 GeV. The result of our numerical computation for the cases 1, 2, 3 and 
4 are illustrated in Figs. 5, 6, 7 and 8, respectively. 

For the sake of comparison let us recall the decay width formula in the one 
Higgs boson case. The decay width in the minimal standard model (MSM) is given 
through the one- loop order by [12] 

1 m? 



1 mi, 19 3V37r Stt^ 
x|l + i^l#(l6-— + 1^)1^- (37) 

This formula gives us r^^^M^^ ^ W^W^) ~ 7.6 GeV for mn = 300 GeV. Note 
that the radiative correction in (37) is as small as 1.3 %. This is due to the smallness 
of the coefficients of mjf/v'^ in (37). 

Let us now look at Figs. 5,6,7 and 8 closely. We can immediately see conspicuous 
differences in the magnitude of computed values in these four cases. The cases 3 and 
4 (Figs. 7 and 8) show that the computed width ranges from 1 to 8 GeV, depending 
on 771(3 and 777^- In the case 1, on the other hand, the width is smaller by one or 
two orders of magnitude or perhaps more depending on rriG, compared to the cases 
3 and 4. The case 2 also gives small decay width for 400 GeV< mo < 700 GeV. 

This situation may be understood in the following way. As mentioned before, in 
the cases 1 and 2, the tree level coupling vanishes and the predicted decay width 
is only due to the radiative corrections. The radiative corrections thus turn out to 
be small for the cases 1 and 2 for the same reason as in the MSM case. This is 
particularlly interesting, if one would notice that the values ma — 900 GeV and 
niA — 1000 GeV are on the verge of the breakdown of perturvative calculation. 

We also note that the 7r7G-dependence of the width for cases 1 and 2 is wild- 
behaved. It changes by one or two orders of magnitude as we vary niG from 400 
GeV to 900 GeV. For large-TTic, the width goes up hke rriQ, which is in accordance 
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with our initial expectation. At any rate, if the width would be this small as in Fig. 
5, the experimental measurement of the width would not be easy. 

Since we have set sin^(Q; — P) — 1 for the cases 3 and 4, the tree-level prediction 
of the width for these cases coincide with the tree-level value in the MSM case 
without radiative correction , which becomes 7.5 GeV. In Figs. 7 and 8, wc sec that 
the radiative correction becomes important for large mo for ttia = 350 GeV and 
ruA = 700 GeV. The corrections reaches 35 % (18 %) for tua = 350 GeV (m^ = 700 
GeV) at ma — 900 GeV. If the uia is as large as 1000 GeV, the radiative corrections 
become comparable with the tree-level calculation for small mo- This is attributed 
to the largeness of niA and shows that the perturbative calculation is barely allowed 
for this choice of itla- The mcr-dependence in Figs. 7 and 8 is not so wild as in 
Figs. 5 and 6. This is due to the fact that the wild ttiq/v^ behavior of the loop 
amplitudes is tamed by larger mc-independent tree amplitudes. 

We have repeated the above calculation by setting ruh — 600 GeV. The results, 
however, do not differ much from those shown in Figs. 5,6,7 and 8. This fact is 
also very interesting, since the value rUh = 600 GeV is close to the boundary of the 
perturbative region. This insensitivity to nih reminds us of the analysis by Bertolini 
[6]. He studied radiative corrections, in the two-doublet model, to the p-parameter 
and to the relation between the Fermi constant and the SU{2)l gauge coupling 
(often denoted by Ar). He argued that, for a = j3, Ar is independent of m/j. 

Suppose that the Higgs boson H has been discovered and the decay width r{H ^ 
Wi^W£) has been measured to a good accuracy. In such a situation, the calculation 
shown in Figs. 7 and 8 may be useful to probe the existence of another Higgs boson, 
i. e. G and A if we know the values of a and p. Of course we do not know at present 
definite value of a or (3, either. We have to look for various means of measuring the 
parameters of the two-doublet model. Being given various measurements, we will 
be able to fit the data in a multi-dimensional parameter space. 
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Finally we add a comment of the dependence of the invariant amplitude 
M.{p^). We have examined M.{p^) numerically in the hope of possibility of future 
confrontation with measurement. With the parameters of cases 3 and 4, this form- 
factor is almost a constant for 200 GeV < < 350 GeV, tuq — 500 GeV, m\ — 
350 GeV, 700 GeV, and 1000 GeV, and docs not show a peculiar behavior. This is 
because the p^-independcnt tree amplitude dominates over the loop amplitudes. In 
the cases 1 and 2, on the other hand, where the tree amplitude vanishes, the form 
factor goes up rapidly as we increase VP^ from 250 GeV to 350 GeV. It magnitude 
is, however, still very small compared to the cases 3 and 4. 



6. Summary and Discussions 

In the present paper, we have investigated the radiative corrections in the two- 
Higgs doublet model from the viewpoint that the loop corrections to the decay rate 
oiH ^ W^W]^ could be a useful probe into the other Higgs bosons. The corrections 
to the w vertex contain terms behaving like rriQ/v^ for large mo- This indicates 
that predicted width could be sensitive to the value ma . We have demonstrated 
the behavior of the width as a function of ma for various choices of mh, mA, oc and 
(5. We have seen to what extent the terms of the power-behavior would become 
important in the decay width formula. 

The calculation presented in this paper can be easily carried over to the decay 
H — > ZlZl, the analysis of which will be given in our future publication [29]. 
The elastic scattering Wi^W£ — > Wi^W£ , ZlZl is also of particular interest. The 
present calculation will become a useful basis for the evaluation of these elastic 
scattering. 

The Higgs boson masses and the radiative corrections in supersymmetric theories 
have been one of important topics in recent literatures [30]. The parameters in the 
Higgs sector in supersymmetric theories are not completely independent, but are 
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constrained by the supersymmetry. Suppose that the decay rate of if — > Wi^W^ 
be evaluated in supersymmetric models. The parameters describing the decay rate 
can not be varied freely. For such a case, it might be that the screening phenomena 
could occur effectively, because of the constrained parameter space. We will come 
to this problem in the future. 
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Appendix A 

Here we summarize our notations for the Feynman integrals corresponding to 
various types of diagrams in Fig. 2, which are expressed by 



fi{m ) = /i' 



(27r)^F-m2 
, 2 ^ , 



d^k i i 



{2n)^ k"^ — m\{p — ky — rn? ' 



(3^ 



Here // is the scale parameter of the D-dimensional regularization. These notations 
will be used extensively in Appendices B,C, and D. If we set m2 = in Eq. (38), 
we obtain 

1 2 

+(i-f)Mi-|;)}. 



(39) 



The formula (38) is also simplified if we put mi = m2 = m, i.e., 

222^ lr2 „ 



h{p\m\m') = ——{—— + 7£- 2 + ln 



I Jl-4mVp2 + l 

V ^l-4mVp2-l 



(40) 



The vertex integral corresponding to Fig. 4 is expressed in terms of the following 
function 

g{p , m„ m„ ^3) = / _ ^^^2 _ ^2 + -mlk^-ml ^^'^ 
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It has been known that this integral is given by a combination of the Spence function; 
g(j)'^,ml,ml,ml) = \ ln(- — 3^){ln(a - ^+) + ln(a - 

+ (4^^-SP(^) + Sp(^)-Sp(^) 
+ Sp(^) + Sp(^)-Sp(^)}. (42) 
Here the Spence function is defined as usual by 

Sp(^) = - / _ln(l-x). (43) 

«/ Ob 

Various quantities appearing in (42) are given by 

m| — + ± (ml — mf + p^)^ — Am^p 

^-^Z^^ "2 - (45) 



^± = ' (44) 



9 ' 2 2 ■ 

mf — 777.3 



Appendix B 

The self-energy of the Nambu-Goldstone boson n^„^(p^), and the two-point func- 
tion Uujoip^) are obtained by evaluating Fig. 2, whose internal particles are specified 
in Table 1. Straightforward calculations show that these functions are expressed as 
in (24) and (25), where we have introduced following functions 

tiUp') = ^cos2(a-/3)/2(p^m^,0) 



4 

+-^sm\a-(3)Mp',mjj,0) 

(777^ — 777q)^ 2/ «\/-/„2 2 2 





2 

V 




-mlf 




y2 




-mlY 


' 2 

V 



sin {a- (5)f2{p ,m^,ma) 
■ cos^(q; - /3)/2(p^ m|, 777^) 
/2(p^m^,mc), 
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(46) 



■sm(2a-2p)f2(p\mjj,0) 



2^2 

m? — TTir . , r 9 / sin^ 3 cos a cos^ /9 sin a , 
+ ^ , ^ sin a - /? {m^ ^— - — + /i- — 

+2m^ cos(a - /?)}/2(p^ m^, m^) 
m?T — rf?r / r 9 / cos^ /3 cos a sin^ /3 sin a , 

sm (3 cos p 



Appendix C 



(47) 



The self-energy of the H boson is decomposed into three parts 

nHHip") = ng), + ngi, + ngg, (48) 

which correspond to Figs. 2(a), 2(b) and 2(c), respectively. These three diagrams 
give us the following results: 

(1) _ 9m|, Wa sin=^a 2. . 2 ^2 ^2 ^ 

'^'"^^ + "^'^^'^^^"'"^^cos^(a-/3)/2(/,mLO 



sin 2/3' 

(2m^ + m|^)%sin2a 2 -2/ 2 2 2^ 

+^sin^(a-/3)/2(/,0,0) 

.m?T ,cos^ /3cosq; sin^/^sino;, 2m^ , ^,,9 

+ — —o a - /3 f 

V sm/5 cosp 

x/2(p^mg,mg) 

1 ,m?. ,cos^/3cosa sin^/^sina, 2m4 , ^,,9 

+ 7^ — ( r-T^ ^-Ti ) ^sin(Q; - d)Y 

2^ V ^ sin (3 cos/3 ^ v ^ ^'^ 

xf2{p^,m\,m\) 
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2 \2 



+ ^±!L_!^ cos^(« - (3)Mp', ml 0), 



(49) 



„C2i ,ml, cosasind sinacos/?. 



f 2 sin /3 cos P sin /? cos /3 



+^sin^(a-/3)]A(m^) 

.m? .cose 1 sin 7 sin ci cos /3, 
+ brl + ■ 2^ sin2acosa-/j 

m?T,cos^Q; sin^ a. .cosacos^ /3 sinasin^/?, 

H — ^( )( 

sin/? cos/3 sin/3 cos/3 ' 

+^sin^(a-/3)]/i(mi) 



rm? ,^,sin2Q;,2 . o, n\ sin2Q;, 
"^ff r^/Sin2Q;, r, ^, sin2Q;,, ^ , 

,3m? ,sin^ acosa sinacos^a.n Sm^r.sin^a cos^o;,,, 

+ ^ ) H ^( ) 

^2f2^ cos/3 sin/3 ' 2^;^^ cos/3 sin/3 

x/i(m^), 

(50) 

/a) 1, cosset sin^ct,^ lsin2Q; . , . 

^Wh = - — ^ ^)Th + --—5 cos a - P)n. 51 

■u sm /3 cos /3 sm 2/3 



Appendix D 

The mixing diagrams between H and h are also of the types of Fig. 2. They are 
again decomposed into three terms 

n,H{p') = n« + nl^l + n^^l , (52) 
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where each term in (52) are given by 



n 



(1) 

hH 



1 



, sin 2a , 



3ml 



cos^ a sin^ a sin 2a . . 
cosp smp sin2p 



3mir , 9 1 9,,cos^q; sin^a,sin2Q; 



ff / 9 J- 9\ / 

Ji{ml + -ml){-: 



sin/3 cos/3 sin 2/3 



)-. — T-^ cos(q; — /3) 



x/2(p^m^,m^) 
^"'^^^sin(2a-2/3)^(p^0,0) 



4^2 

, m? , sin^ /3 cos a sin a cos^ /3 , 2m^ 

+{— + + 

V cos p smp V 

2ml 



cos q; 



m% . cos^ /? cos a sin a sin /3 , 



sin /3 



cos /3 



sin(Q; — /3)} 



x/2(p^m^,m^) 

1 , m? , sin^ /3 cos a sinacos^/^, 2m 4 , ^, , 

2 V cos/3 sin/3 v 



, cos^ /3 cos a 



V sin /3 
xf2(p'^,m\,m\) 
1 

' 



sinQ;sin^/3, 2m 4 . , 

^ ^sm a-/3 } 

cos /3 v 



-m2.)sin(2a-2/3)/2(p^m2,,0) 



1 



2t;2 



(m^ - m\){m]j - m\) sm{2a - 2/3)/2(p^ m\, 0), 



(53) 



n 



(2) 
hH 



m\ cos a sin /3 sin a cos /3 



-2t;2' 



cos^/3 



+ 



sin2/3 



-) sin 2q! sin(Q! — /3) 



m|^ cos a cos /3 sin o; sin (3 



2v' 



-( 



sin /3 



cos^ /3 



) sin2Q;cos(Q; — jS) 



-^sin(2a-2/3)}/i(m^) 



m| , cos a sin /3 sin a cos /3 , . 



+ 



4i>2^ cos^/3 ' sin^/3 
m|j^ , cos a cos /3 sin a sin /3 , . 



) sin2Q;sin(Q; 



Av'^ ^ sin^ /3 



cos^ /3 



) sin2Q;cos(Q; 



/5) 
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-|;|sin(2a-2/3)}/i(mi) 



3m; sin 2a, 
2>m'\j sin 2a cos^ a 



'sin2(Q; - {3) 



sin"' ot . 



2f2 sin 2/3 sin/? cos/9 



+{ 



3m|^ sin2Q;, 
4^2 *^ sin 2/3' 



'sin2(Q; - /3) 



+ 



3m? sin 2q; , cos" 



a 



sin^ a 



cos(a - f3)}h{m]j) 



2v^ sin 2/3 cos/3 sin/3 



)sin(Q;-/3)}/i(m^), 



1 sin 2a 
V sin 2/3 



{cos(q; — P)Th + sm{a — P)Th}. 
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Table 1 



Combinations of internal particles {X, Y) running in Fig. 2(a) and X in Fig. 2 
for m^land U% {i = l,2). 



Propagator 


Internal particle species 


n(i) 

WW 

n(2) 

WW 

tt(1) 


ih,w), iH,w), {h,G), {H,G), {A,G) 
G, A, h, H 

{h,w),{H,w),{h,G),{H,G) 
G, A, h, H 



Table 2 

Combinations of internal particles (X, Y) running in Fig. 2(a) and X in Fig. 2 
for n«^and n^i(^ = l,2). 



Propagator 


Internal particle species 


n(i) n^') 
n(2) n^^^ 


{H,H), {h,H), {h,h), iw,w), {z,z), 
{G,G), {A, A), {A,z), {G,w) 
G, A, h, H 
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Table 3 



Combinations of internal particles (X, Y, ; Z) in Fig. 4(a) and (X, Y) in Figs. 4(b) 
and 4(c)for the vertices F*^*^ (i — 1, 2, 3). 



Vertex 


Internal particle species 


r(i) 


{w,w-h),{w,w-H),{G,G-h),{G,G-H), 




{G,G-A), {A,A-G), (G,w;h), (w,G;h), 




{G,w;H), {w,G;H), {H,H;w), {H,H;G), 




{H,h;w), {h,H;w) {H,h;G), {h,H;G), 




{h, h; w), {h, h] G) 


r(2) 


{w,w), {z,z), {G,G), {A, A), {A,z), 




{G,w),{H,H),{H,h),{h,h) 


r(3) 


{w,h),{w,H),{G,h),{G,H),{G,A) 
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Figure Captions 



Fig. 1 

Tadpole diagrams of H and h fields contributing to (a) Th and (b) T^. 
Fig. 2 

Self- Energy diagrams contributing to (a) 11 -j^ (b) n-^-*, and (c) nj|\ respec- 
tively. The diagram (c) denotes the counter terms, which are expressed by the 
tadpole contributions. A pair of indices (i, j) refers to either of {w, w), {w, G), 
{H, H), or [h, H). Internal particle species are given in Tables 1 and 2. 

Fig. 3 

General configuration of the Hww^ vertex. The Nambu-Goldstone bosons are 
put on the mass shell {p\ = p\ = 0), while the the Higgs boson H are kept 
off-shell to keep generality. 

Fig. 4 

Radiative corrections to the w vertex contributing to (a) F'^^^ (b) F^^^ 
and (c) F^^). Internal particle species {X,Y]Z) in (a) and {X,Y) in (b) and 
those in (c) are given in Table 3. 

Fig. 5 

The decay width (35) as a function of mc- The mixing angles are determined 
by tana = tan/3 = 2. The masses of the neutral Higgs bosons are assumed to 
be niH — 300 GeV, = 400 GeV. The CP-odd Higgs boson mass is taken 
as (a) niA — 350 GeV (sohd hue), (b) = 700 GeV (dashed line) and (c) 
niA — 1000 GeV (dotted hue), respectively. 
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Fig. 6 

The decay width (35) as a function of ma- The mixing angles are determined 
by tana = tan/3 = 10. The masses of the neutral Higgs bosons are assumed 
to be niH = 300 GeV, rrih = 400 GeV. The CP-odd Higgs boson mass is taken 
as (a) niA — 350 GeV (sohd hne), (b) ttia — 700 GeV (dashed line) and (c) 
niA — 1000 GeV (dotted hne), respectively. 

Fig. 7 

The decay width (35) as a function of ttig- The mixing angles are determined 
by tan/? = 2 and sin^(Q; — P) = 1. The masses of the neutral Higgs bosons 
are assumed to be nin = 300 GeV, irih = 400 GeV. The CP-odd Higgs boson 
mass is taken as (a) = 350 GeV (solid line), (b) — 700 GeV (dashed 
hne) and (c) = 1000 GeV (dotted line), respectively. 

Fig. 8 

The decay width (35) as a function of ttig- The mixing angles are determined 
by tan/3 = 10 and sin^(a — /3) = 1. The masses of the neutral Higgs bosons 
are assumed to be uih = 300 GeV, rrih = 400 GeV. The CP- odd Higgs boson 
mass is taken as (a) ttia — 350 GeV (solid line), (b) uia — 700 GeV (dashed 
line) and (c) niA — 1000 GeV (dotted line), respectively. 
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